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Abstract 

We present recent results for the Landau-gauge gluon and ghost propagators in SU (3) lattice gluodynamics obtained on 
a sequence of lattices with linear extension ranging from L = 64 to L = 96 at (3 — 5.70, thus reaching "deep infrared" 
momenta down to 75 MeV. Our gauge-fixing procedure essentially uses a simulated annealing technique which allows us 
to reach gauge-functional values closer to the global maxima than standard approaches do. Our results are consistent 
with the so-called decoupling solutions found for Dyson-Schwinger and functional renormalization group equations. 
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1. Introduction 

The infrared behaviour of gauge-variant Green's func- 
tions of Yang-Mills theories has become an increasingly 
interesting topic during the last decade. The interest was 
originally stimulated by confinement scenarios proposed 
a long time ago by Kugo and Ojima 0, E|, Gribov [|| 
and Zwanziger 0, Q. The recent progress is due to the 
discovery of consistent asymptotic solutions of the whole 
tower of Dyson-Schwinger (DS) equations and more re- 
cently of functional renormalization group (FRG) equa- 
tions in the deep infrared (IR) limit [& 0, S, H, E3, LLL, El | . 
These solutions called scaling or conformal solutions be- 
have accordin g to power laws with well-determined ex- 
ponents [1, [9l. Txxl| - The expectation is that they respect 
global BRST invariance. Since these solutions lead to a 
vanishing gluon propagator and correspondingly to an IR- 
singular ghost dressing function they fit nicely with the 
aforementioned scenarios. The running coupling related to 
the ghost-ghost-gluon vertex exhibits an IR fixed point Q 
as also proposed by Shirkov and Solovtsov [lJ [l4| ■ 

There is a different set of so-called decoupling solutions 
as proposed in [H, EE 12, IS. flit- These solutions - 



well as providing the expected Polyakov loop behaviour 
at the deconfinement transition in pure SU(2) and SU{2>) 
gauge theories [2l]]. Moreover, both types of solutions, 
when interpolated from the infrared asymptotics to the 
perturbative region by solving numerically the (properly 
truncated) system of DS or FRG equations, behave quite 
similarly in the momentum range relevant for hadron phe- 
nomenology. 

Even if it might appear to be an academic question, it 
is important to ask which set of solutions is the correct one. 
Ab-initio lattice gauge theory computations are expected 
to solve this issue. SU(2) and SU(3) lattice computations 
of Landau-gauge gluon and ghost propagators have been 
carried out by several groups. For a recent review see 
22j and papers cited therein. On the one hand, as long 



12J — are characterised by 577(3) see [13, [24 , 



thoroughly discussed also in 

a non-zero IR gluon propagator as well as by an IR-finite 
ghost dressing function, i.e., they do not agree with the 
Kugo-Ojima criterion. The name decoupling reflects the 
fact that the corresponding running coupling decreases to- 
wards zero in the limit of vanishing momenta. But this 
behaviour does not mean that the decoupling solutions 
contradict gluo n and quark confinement (see also the dis- 
cussion in [20|). Both sets of solutions demonstrate the 
expected positivity violation of the gluon propagator as 



as in the four-dimensional case the linear lattice sizes did 
not reach far beyond 0(5 fm) the ghost dressing function 
was observed to rise towards the infrared limit. When 
fitted with a power law the corresponding exponent was 
found to be much smaller than predicted by the scaling 
solutions. On the other hand, the gluon propagator has 
not yet been found to decrease towards the IR limit but 
rather approaching a finite plateau value at p — (for 



25(). However, investigations of the DS 
equations on a 4D torus have demonstrated that linear 
box sizes of O(10 fm) or even more, might be necessary in 
order to see the correct asymptotic behaviour. 

In the meantime simulations of SU(2) pure lattice gauge 
theory have reached huge lattice sizes [22, [27| confirming 
nothing but an IR-plateau behaviour for the gluon prop- 
agator. Moreover, in [28j Cucchieri and Mendes derived 
infrared bounds for the gluon and ghost propagators indi- 
cating, e.g., a non-zero value for the zero-momentum gluon 
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propagator in the infinite-volume extrapolation. 

However, the aforementioned results of DS and FRG 
studies as well as most of the lattice results obtained with 
the standard minimal Landau gauge did not take into 
account the effect of Gribov copies. For SU(2) some of 



us have shown 29, 30, 31, 32, 33] that the influence of 



Gribov copies is more severe than many authors might 
have expected (see also the recent investigation in [34|). 
The effect is even stronger if non-periodic Z{2) gauge- 
transformations (flips of all link variables orthogonal to 
fixed 3D sheets with a factor -1) are taken into account. 
Already on modest lattice volumes the gluon propagator 
as well as the ghost dressing function were seen to run into 
IR plateaus indicating that lattice results seem to s upp ort 
the decoupling solution of DS and FRG equations [331 ]. 

In the SU(3) case the effect of Gribov copies on the 
Landau gauge gluon and ghost propagators has been stud- 
ied without Z(3) flips so far [23|, |35j, [36J showing that only 
the ghost propagator seems to be systematically affected 
within 5 to 10%. 

In the meantime, we have extended the computation 
of Landau gauge gluon and ghost propagators in SU(3) 
gluodynamics to linear lattice sizes of around 0(16 fm). 
Preliminary results for lattice sizes up to 80 4 have been 
presented already in [37| clearly showing that an infrared 
plateau of the gluon propagator evolves and demonstrat- 
ing — to our knowledge for the first time — a flattening 
of the ghost dressing function as well. Here we go a step 
further in order to get more confidence by increasing the 
lattice size for the gluon propagator up to as much as 96 4 
as well as by increasing considerably our statistics for both 
propagators. In our investigation we have put some em- 
phasis on careful gauge fixing with the simulated anneal- 
ing method which wins in efficiency in comparison with 
standard over-relaxation, the more degrees of freedom the 
system has (see below and discussions in [3l|). 

2. General setup 

We compute the SU(3) gluon and ghost propagators 
with Monte Carlo (MC) techniques on a lattice with pe- 
riodic boundary conditions. The standard Wilson single- 
plaquette action and the lattice definition for the gauge 
potentials 



A^x + fi/2) 
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are adopted. In order to fix the Landau gauge for each 
lattice gauge field {U} generated by means of the MC pro- 
cedure, the gauge functional 



! 4 

p u\9] = 5 E E ^ e Tr 9* u *» gt +ft 



(2) 
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is iteratively maximised with respect to a gauge trans- 
formation {g x } which is taken as a periodic field as well. 
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Figure 1: Example of a SA schedule. Shown is the temperature 
step size AT by which T is reduced in the SA algorithm used on a 
96 4 lattice (see text for further details). 



In order to approach the global maximum (related to the 
fundamental modular region) as closely as possible, we use 
the simulated annealing (SA) algorithm 38], in combina- 
tion with subsequent standard over-relaxation (OR). The 
latter is applied in the final stage of the gauge-fixing proce- 
dure in order to finalise the transformation to any required 
precision of the transversality condition V^yLj = 0. More 
than a decade ago, SA was shown to be very efficient when 
dealing with the maximally Abelian gauge 3(1, 4(3] ■ 

In case of the present approach the SA algorithm gen- 
erates a field of gauge transformations g x by MC iterations 
with a statistical weight proportional to exp(Fu[g]/T). 
The temperature T is an auxiliary parameter which is grad- 
ually decreased in order to maximise the gauge functional 
Fu[g]- In the beginning, the T chosen must be large 
enough to allow the traversing of the configuration space 
of g x fields in large steps. An initial value T init = 0.45 



was found to be sufficient for that [4l|. After each sweep, 
including one heatbath and four micro-canonical update 
steps at each lattice site, T is decreased until g x is cap- 
tured in a particular basin of attraction. We choose the 
lowest temperature value to be Tfi na i = 0.01 and use a fine- 
tuned simulated annealing schedule before applying final 
over-relaxation steps to reach Landau gauge with good 
precision (max. JHe Tr^A^V^^)] < 10~ 13 ). 

An infinitely slow and long simulated annealing process 
would definitely lead to the global maximum. In prac- 
tice, however, we chose schedules according to which the 
temperature is reduced in steps of finite but varying size. 
These schedules are motivated by the following observa- 
tion. When the thermal average (Fjj[g\) is considered as 
a function of T, one easily notices its monotonous but 
non-linear dependence. In fact, (Fjj[g]) increases when 
T decreases, but its gradient reaches a strong maximum 
at a "critical" value T ~ 0.405 ± 0.01, which resembles 
very much a phase transition 4l|. This means that after 
starting with small AT steps, as required for proper thcr- 



2 



malisation close to Tj n it = 0.45, it is important to keep 
the step size especially small within the narrow critical re- 
gion T ~ 0.40 . . . 0.41. Indeed, the distribution of the final 
gauge-functional values (those of the local maxima after 
simulated annealing and over-relaxation) is shifted to no- 
ticeably higher values if the temperature step size AT is 
reduced according to the rise of (Fu[g]), i.e., if AT is taken 
relatively small within the critical interval. This is in par- 
ticular an efficient approach, as almost no improvement 
was observed, when further reducing AT outside this in- 



terval [41 1 



As an example, Fig. [T] shows the step sizes we use for 
the different temperature intervals on a 96 4 lattice, where, 
especially around T ~ 0.40... 0.41, AT is taken to be 
very small. Such extremely long SA "chains" , with an or- 
der O(10 4 ) of iterations, allow us with high probability to 
reach local maxima of the gauge functional Fjj [g] close to 
the global maximum (as far as possible for the computer 
resources available), i.e., to the fundamental modular re- 
gion, even with only one gauge-fixing attempt ("first-copy 
approach"). 

Note that we do not apply here Z(3) flips, the so-called 
FSA method [30j, |32| which in principle is capable of pro- 
viding even larger Fjj[g\ values than our SA algorithm 
does. But we have seen the difference between SA and 
FSA results decreases when L increases [Hj]. 

The computations presented are carried out at rather 
strong coupling, /3 = 6/<?q = 5.70. The reason for this 
choice is to get access to comparatively large physical vol- 
umes at the price of a rather coarse lattice (a « 0.17 fm). 
In order to study the volume dependence we calculate the 
gluon propagator for linear lattice sizes L ranging from 
64 to 96. Thus, our largest lattice size corresponds to 
(16 fm) 4 . The ghost propagator is studied for linear lat- 
tice sizes L — 64 and L = 80. At larger lattice sizes and 
the lowest momenta the inversion of the Faddeev-Popov 
matrix turns out to converge only in rare cases. Obvi- 
ously, this is due to the occurrence of very small eigen- 
values which generate some algorithmic problems. There, 
a modification of the used matrix inversion method (see 
Sec. or an even better gauge fixing, driving configura- 
tions further away from the Gribov horizon [36j , could be 
valuable means of reducing such problems in future, but 
it is beyond the scope of the current work. 

3. Gluon propagator 

The gluon propagator is defined by 



DfM = (A«(k)AU-k)) 

= (s, v -^)s ab D(q 2 ). 
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Figure 2: The bare lattice gluon propagator D(q 2 ) versus q 2 for 
/3 = 5.70 and various lattice sizes. We also show data on D(0) (left). 

with fc M e (-L/2, L/2}. For q ^ 0, one gets 



8 4 



(4) 
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whereas at q = the zero-momentum propagator D(0) is 
defined as 



(5) 



a=l fj,= l 



Our data for the gluon propagator obtained for various 
lattice sizes is presented.!] in Fig. [H One can clearly see 
a flattening of the gluon propagator as a function of q 2 
for small momenta. Note also the weak volume depen- 
dence of the results. To illustrate the latter, in Fig. [3] we 
present also the dependence of the zero-momentum prop- 
agator D(0) (acc. to Eq. $5§) on the inverse lattice size 
1/L. From this point of view it is natural to conclude 
that in the infinite volume limit the gluon propagator will 
approach some constant value as q 2 — > 0, i.e., the gluon 
dressing function Z(q 2 ) = q 2 D(q 2 ) seems to decrease lin- 
early with q 2 . The observed IR behaviour related to the 
zero- momentum modes of the lattice gauge field (x) can 
be associated with a massive gluon. 

Our SU (3) IR gluon plateau results are in close agree- 
ment with analogous SU(2) results found recently on huge 
AD symmetric lattices [12, Hi]- 



(3) 



1 In this letter the gluon and ghost propagator data has not been 
renormalised in contrast to our former studies, in particular in |37l| . 



where A(k) represents the Fourier transform of the gauge 
potentials according to Eq. ([T]) with Landau-gauge fixed 
links. The momentum q is given by q^ — (2/a) sin (nk^/L) 
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Figure 3: Zero-momentum gluon propagator D(0) versus 1/L. 
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Figure 4: Bare ghost dressing function J(q 2 ) versus q 2 for L = 64, 80 
at (3 = 5.70. Errors are not shown at the two lowest q 2 (squares). 



4. Ghost propagator 

The Landau-gauge ghost propagator is defined by 
G a \q) = a 2 ^2(e-^ lk < x -y^ L [M- 1 ]f v 



<*° W) , 



(6) 



where M denotes the lattice Faddeev-Popov operator, be- 
ing the Hessian of the gauge functional ^ with respect to 
g x , in the background of the gauge- fixed links U xp , 



with 
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T a , a — 1,...,8 are the (Hcrmitian) generators of the 
su(3) Lie algebra satisfying Tr [T a T b ] = 6 ab /2. 

To invert M we use the conjugate gradient (CG) al- 
gorithm with plane-wave sources ?/> c with colour and po- 
sition components ip^(x) — 5 ac exp(2ni k-x/ L). In fact, 
we apply a pre-conditioned CG algorithm (PCG) to solve 
M£ b (f> b (y) = xp2(x) where as pre-conditioning matrix we 
use the inverse Laplacian operator A -1 with diagonal colour 
substructure (for details see 23, 42|). 

For the large lattice sizes as considered here, we are 
confident that finite- volume distortions for all lattice mo- 
menta besides the two minimal ones do not change consid- 
erably with increasing L (see Fig. 0] and [13] for details). 
In this figure the ghost dressing function J(q 2 ) = q 2 G(q 2 ) 
is presented in a log-log scale. We do not see any power- 
like singular behaviour in the limit q 2 — > 0. Instead, we 
have a good indication that J(q 2 ) reaches a plateau just 
as the de coup ling solution of DS and FRG equations does 
(see also [ml^y 



5. Running coupling 

Finally, let us present the running coupling defined as 
the renormalization group (RG) invariant product 



a s (q 2 ) = ^Z(q 2 )J 2 (q 2 ), 



(8) 



of the Landau-gauge gluon and ghost dressing functions. 
This definition is based on the ghost-gluon vertex in a 
momentum-subtraction scheme with the vertex renormal- 
isation constant (in Landau gauge) set to one. This is 
possible [43j], since the vertex is known to be regular 44 1 



(see also the lattice studies [4a, |46[). Note that the rela- 
tion of a s in this scheme to the running coupling in the MS 
scheme is known to four loops and it can provide a valu- 
able alternative to the MS coupling in phenomenological 
applications [43 1. 

Beyond perturbation theory, the behaviour of a s dif- 
fers at low scales for the scaling and decoupling solutions. 
Based on our propagator data we can calculate a s for in- 
termediate and lower scales, and it clearly shows a decrease 
towards q 2 — > (see Fig. [5]). This is again consistent with 
the decoupling DS and FRG solutions. 

6. Conclusions 

The progress achieved on the lattice during the last 
two years in studying the IR limit of gluodynamics and 
checking the well-known scenarios of confinement in terms 
of Landau-gauge Green's functions leads us to the follow- 
ing conclusions. Within the standard lattice approach as 
described above only the decoupling-type solution of DS 
and FRG equations seems to survive. Since for this solu- 
tion the gluon propagator tends to a non-zero IR value, it 
corresponds to a massive gluon. It has been argued that 
this behaviour contradicts global BRST invariance [l2j ]. 



4 



9 

Co 

3 




100 



Figure 5: Running coupling a a (q 2 ) versus q 2 for lattice sizes 64 4 and 
80 4 at = 5.70. 



But the lattice approach as discussed here has a few 
weak points. The choice of the gauge potentials A^{x) 
and correspondingly of the gauge functional Fjj(g) is in no 
way unique. As long as we are reaching the infrared limit 
by employing rather large lattice couplings the continuum 
limit is not under control. Moreover, we have used stan- 
dard periodic boundary conditions which certainly have 
an impact on the IR limit. The fact that under these 
conditions the gluon propagator does not tend to zero is 
related to the behaviour of the zero-momentum modes, 
which do not become sufficiently suppressed as the lattice 
size increases. Changing the definition of lattice Landau 
gauge, and correspondingly the lattice definitions of (x) 
and M, modifying the boundary conditions and further 
improving the gauge-fixing procedure, e.g., by taking Z(3) 
flips as mentioned in Sec. [1] into account, may essentially 
suppress the zero-momentum modes and correspondingly 
change the behaviour of both the gluon and ghost propa- 
gators. Therefore, a final conclusion still cannot be drawn. 

Note that one of us (A.S.) has recently carried out a 
lattice computation in the strong-coupling limit. For the 
gluon and the ghost propagator at larger a 2 q 2 it was pos- 
sible to extract the right exponents as expected for the 
scaling solution [47| . At asymptotically small momenta, 
however, results were shown to depend strongly on the 
lattice definitions of A^{x) and M. This corresponds to 
observations in studies of DS and FRG equations, namely 
that for decoupling-like solutions any IR-asymptotic val- 
ues of the gluon propagator and the ghost dressing function 
can be considered as boundary conditions [lij ]. 



It has been argued in 



that a BRST-invariant gauge- 
fixing prescription is possible on the lattice. It remains 
to be seen, whether the preferred scaling behaviour of 
Landau-gauge gluon and ghost propagators can be achieved 
consistently also in lattice Yang-Mills theories. 
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